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Abstract 

Predictions for the angular distribution of massive quarks and leptons are presented, including 
QCD and QED corrections. Recent results for the fermionic part of the two-loop corrections to the 
electromagnetic form factors are combined with the BLM scale fixing prescription. Two distinctly 
different scales arise as arguments of a s (fi 2 ) near threshold: the relative momentum of the quarks 
governing the soft gluon exchange responsible for the Coulomb potential, and a large momentum 
scale approximately equal to twice the quark mass for the corrections induced by transverse glu- 
ons. Numerical predictions for charmed, bottom, and top quarks are given. One obtains a direct 
determination of av(Q 2 ), the coupling in the heavy quark potential, which can be compared with 
lattice gauge theory predictions. The corresponding QED results for r pair production allow for a 
measurement of the magnetic moment of the r and could be tested at a future r-charm factory. 
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1. Introduction 



QCD predictions for the angular distributions of partons can be tested in a variety of ways. For 
massless partons, i.e. light quarks or gluons, the canonical approach is based on angular distribu- 
tions of jets or related quantities, like the thrust axis. In fact, the observation of a distribution 
cx (1 + cos 2 6) for two-jet-events has provided convincing evidence for the spin- 1/2 nature of partons 
produced in e + e~ annihilation. Alternatively, for charmed or bottom quarks one may directly iden- 
tify the heavy quark and heavy meson directions. This approach has been employed frequently in 
the analysis of the b quark angular distributions at high energies. In this kinematical region gluon 
radiation affects the shape of the distribution, typically rendering the lowest order predictions more 
isotropic. Real and virtual corrections must be considered jointly, since only the sum leads to an 
infrared finite prediction. 

In this paper a different kinematical region will be considered: heavy quark production close to 
threshold, where the cross section is dominated by a few exclusive channels. In this case one expects 
that the heavy mesons will essentially follow the heavy quark direction — a simple consequence of 
Newton's law of inertia. To employ the parton model in the case of hadron pair production may be 
somewhat surprising, since the structure of form factors for exclusive channels imposes stringent 
constraints on the angular distributions of the hadrons. However, as a consequence of the near 
degeneracy of the pseudoscalar and vector mesons B and B* , the sum of the various channels (BB, 
BB + c.c, B*B ) may easily combine to saturate the prediction for the angular distribution of the 
quarks. For definiteness, we will subsequently discuss the case of b quarks. For certain kinematical 
regions our calculation can be applied for charm quarks and will be valid a forteriori for top quarks. 

Close to threshold, in the limit fi = \q \/ ^Jq 2 + vOq = ^Jl — 4toq / s —> 0, the cm angular distri- 
bution for e + e~ — ► Q Q is of course isotropic, a result of 5-wave dominance. The small admixture 
of P- waves slightly above threshold provides a contribution oc /3 2 cos 2 6 which will be studied in 
this paper. The tree graph prediction for the angular distribution is trivial. The leading 0(a s ) 
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corrections have been evaluated in [1] for arbitrary quark mass and energy. However, unavoidably, 
the renormalization scale at which the strong coupling a s (ji 2 ) has to be evaluated in these formulae 
can only be fixed by a two-loop calculation. In the high energy region one may convincingly argue 
on physical grounds that the scale for the dominant correction terms should be of order y/s. In 
the threshold region, however, both y/s and the quark momentum \q\ are viable options for the 
energy scale. In the BLM procedure [2] the renormalization scale is fixed by resumming all terms 
proportional to the QCD beta-function into the running coupling. To carry out the BLM procedure 
to leading order, it is sufficient to identify the raj terms in the next to leading order coefficients. 
Recent results for the fermionic part of the two-loop corrections to the total rate [3], and in partic- 
ular to the 7* — ► Q Q form factors allow for a solution of the scale-setting problem, at least in the 
threshold region. In fact, employing the BLM scale-setting prescription, one may even fix not only 
the scale but perhaps even accommodate the bulk of the two-loop corrections. 

Real radiation is strongly damped in the threshold region, decreasing oc (3 2 relative to the tree 
contributions. Hence, dropping the terms of higher power in /3, one may dispose of real radiation 
(even if multiplied by an infrared cutoff) and limit the discussion to virtual radiation; i.e. to the 
corrections to the form factors. Contributions from the instantaneous Coulomb potential and from 
hard photon exchange can be clearly distinguished, allowing for the resummation of the former. 
The contribution from the Coulomb part of the heavy quark rescattering yields a series in (ira s / (3) n . 
If one utilizes the ay scheme, defined for the heavy quark potential V{Q 2 ) = —4ir Cp ay(Q 2 )/Q 2 , 
then we find that the scale of ay is set to a value proportional to s /3 2 , the square of the relative 
momentum. Thus one can obtain a determination of ay(Q 2 ) from measurements of e + e~ —^QQ 
near threshold. In particular we shall show that the "anisotropy" A(p 2 ) defined by 

-ocf + Acos 2 6> (I) 

d cos 

is sensitive to of ay(s (3 2 ). Precise predictions for ay(Q 2 ) have been given in [4] using heavy quark 
lattice gauge theory with constraints from the T spectrum. 

We also will discuss the application of the anisotropy formalism to the QED process e + e~ — ► 
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t + t~ . We shall show that the anisotropy of the angular distribution provides a new way to measure 
the Pauli form factors of the leptons, the analog of the anomalous magnetic moment in the timelike 
region. In addition, the radiative corrections we discuss here change the normalization of the cross 
section near threshold. 



2. Form factors, angular distributions, and infrared singularities 



The amplitude for the creation of a massive fermion pair from a virtual photon is characterized by 
the Dirac (F\) and Pauli (F 2 ) form factors: 

uA^v = ieQfuij^F^q 2 ) + q" F 2 (q 2 )]v (2) 

where a^ v = ^[7^,7;/]- The photon momentum flowing into the vertex is denoted by q, the fermion 
mass by to. The resulting angular distribution is conveniently expressed in terms of the electric 
and magnetic form factors G e and G m [5]: 



dcr(e+e- -+ff) _ a 2 Q}fi 



dfi 



4s 



4 to 



\G e \ 2 sm 2 + |GU 2 (l + cos 2 0) 



with 



G e = F 1 

The anisotropy is thus given by 

A 



4 to' 



;F 2 , 



G r 



F 1 + F 2 . 



(4) 



\G r 



(l-/3 2 )|G e 



IC^I 2 + (1 - /3 2 )|C e | 2 
A 



1 - A 



where 



A = — 



fi 2 \ Fl \ 2 (l-fi 2 )-\F 2 \ 2 



(5) 



(6) 



2 \F X + F 2 \ 2 (1 -/3 2 , 

In Born approximation Asom = /3 2 /(2 — (i 2 ). Note that for F 2 = 0, the anisotropy is identical 
to the Born prediction, independent of F\. Thus the form 

A 



1 



■2F,(s) 



1+o 

7T 



A -t 

^iBorn — 2 
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isolates i*2( s )- This provides a way to experimentally determine the timelike Pauli form factor of 
the t lepton. The QED prediction is 

which is, up to the sign, equal to the familiar Schwinger result ^(0) = a/2ir. Away from threshold 
the one-loop QED prediction is 

1-/3 2 



F 2 ({3) = - 

7T 



■tn- 



(9) 



4/3 1 + /3 

We neglect this type of higher twist corrections in the following. 

In Born approximation F\ = 1 and _F 2 = 0. The impact of one- and two-loop radiative cor- 
rections on the form factors and angular distributions will be discussed in the context of QED in 
the remainder of this section; the case of immediate interest, namely QCD, will be discussed in 
section 3. To demonstrate the line of reasoning, we first shall present the arguments for the leading 
order calculation in some detail. 

In the QED calculation the order a correction to the Dirac form factor 8F\ in the timelike 
region exhibits an infrared singularity which can be regulated using a nonvanishing photon mass A: 

8F r = 6F? n + 6F* R tn^ (10) 

with 

«f = f£-!f+^ + <^», (id 

«f = -\-f- + 0(l3 l ). (12) 
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The leading term of is proportional ira/ fi and exhibits the familiar Coulomb singularity. Also 
the constant term and the term linear in /3 are infrared finite. The infrared singular part of F\ is 
strongly suppressed at threshold oc /3 2 , giving rise to a /3 3 contribution to the rate. The correction 
to the Pauli form factor 8F2 is infrared finite and approaches a constant value at threshold: 

6F 2 = -\- + 0(p 2 ). (13) 

I 7T 



Real radiation, in contrast, vanishes as (3 3 in the threshold region, where two powers of (3 result 
from the square of the dipole matrix element, and one power of (3 comes from phase space. It 
exhibits the same logarithmic dependence on the infrared cutoff as the F\ form factor and the same 
leading (3 dependence as the infrared singular part of the virtual correction. As a consequence 
of the strong suppression oc (3 3 it can be neglected in the threshold region, together with the 
corresponding infrared divergent part of the form factor. The angular distribution and, similarly, 
the correction to the total cross section in the threshold region are therefore determined by the the 
infrared finite parts of the form factors. To order a one thus finds for the coefficient describing the 
angular dependent piece 

■ 4 =-M i+ 4^)- (14) 

As we shall show there are interesting modifications of the anisotropy due to the running of the 
QCD coupling, and the dependence of the renormalization scale on y/s and \q \ will be crucial. 

The (9(a 2 )-QED corrections to the form factors, induced by light fermion loops, have been 
calculated analytically in [3]. In the threshold region one obtains 




The calculation has been performed in the limit where the mass of the light virtual fermion mj is 
far smaller than to, a situation appropriate for the subsequent translation to QCD. The factor raj 
is introduced to allow for several light fermions and, in our case, to single out the fermion-induced 
terms. These formulae provide the first step on the way to a full two-loop calculation in order 
a 2 . As we shall see, the results require two conceptionally different scales in the argument of the 
running coupling, a scale of order s from the hard virtual correction from transverse photons and 
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a soft scale of order s /3 2 from the Coulomb rescattering. Supplemented by the BLM prescription 
they even determine the dominant two-loop gluon-induced terms in QCD. 

The linear combination appearing in the denominator of A in Eq. (6) is thus given by 



1 



a ir 
4/3 

a 



1 



1 



i 



V- [in- , 



8 = 1 



i) 

11 N 

12 



(17) 



The rif terms arise from the vacuum polarization insertions and thus can be resummed into the 
QED running coupling: 



a(Q 2 



a 



1 



:i8i 



The constant 5/3 is the usual term in the Serber-Uehling vacuum polarization II(Q 2 ) at large Q 2 . 
This corresponds to the usual QED scheme where V{Q 2 ) = — 4ir a(Q 2 ) / Q 2 is the QED potential 
for the scattering of heavy test charges. One thus obtains 



1 «(^ 2 )vr 2 a(se^lA) 



1 



4/3 

a(se 3 / 4 /4)' 



1 



a(s (3 2 ) 7r 
4/3 , 



(19) 



Two distinctly different correction factors arise. The first originates from hard transverse photon 
exchange, with the scale set by the short distance process; the second from the instantaneous 
Coulomb potential. It is remarkable and non-trivial that the non-logarithmic terms in the ira/(3 
corrections are absorbed if the relative momentum is adopted as the scale for the coupling. Up to 
two loops the running coupling governing the Coulomb singularity is thus identical to the running 
coupling in the potential. This will provide an important guide for the application of these results 
to QCD. 

The proper resummation of the 1/(3 terms based on Sommerfeld's rescattering formula then 
leads to 



l^i + F 2 



1 



a(TO 2 e 3 / 4 )\ x 



1 - e- 



f20) 



with 



x 



a(4m 2 /3 2 ) ir 

is 



(21) 



In a similar way one finds for the relevant combination in the numerator of (6) 



F t \ 2 - \F 2 \ 2 ^ 1 - 3 



a 




IT 



) 



) 



X 



(22) 



with 



x 



a(4m 2 (3 2 /e) ir 

is 



(23) 



The I-F2I 2 term in the numerator can actually be ignored in the present approximation. The scales 
of the effective coupling differ between the numerator and denominator of (6): In particular in the 
factor arising from Coulomb exchange the scale is significantly smaller in the numerator than in 
the denominator. This behavior is consistent with qualitative considerations based on the relative 
distances relevant for S- versus P- waves in the Coulomb part. In the factor arising from hard 
photon exchange the scales are quite comparable, with a slightly larger value in the numerator. 

One thus arrives at the prediction in the context of QED for the anisotropy which involves four 
scales: 



To display the effects more clearly, the ratio of the anisotropy to the Born prediction A/A^ OTn is 
shown in Fig. 1 for the case of r pair production. The dashed curve gives the prediction for constant 
ciqed; the solid curve shows the effect of the lepton vacuum polarization Tl(Q 2 ) in the QED running 
coupling. The vacuum polarization affects the anisotropy for small (3 because two different scales 
appear in the S- and P-wave Coulomb rescattering corrections. Away from threshold A essentially 
measures the anomalous magnetic moment. 

3. From QED to QCD 

The QED coupling a(Q 2 ) translates into the the QCD coupling ay(Q 2 ), defined as the effective 



A = 



1-1' 



A 




(24) 
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A/A 



Born 
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Figure 1: Ratio between the anisotropy A and the Born prediction Asom as function of (3 for the 
process e + e~ — ► r + r~. Dashed curve: constant a; solid curve: including the running of a. 



charge in the potential 



4 4Tra v (Q 2 ) 
3 



(25) 



for the scattering of two heavy quarks in a color-singlet state. In the BLM procedure all terms 
arising from the non-zero beta-function are resummed into ay(Q 2 ). For example, all raj-dependent 
coefficients vanish in the ira/(3 terms if the scale of the relative momentum is adopted. This is, 
in fact, a result expected on general grounds: threshold physics is governed by the nonrelativistic 
instantaneous potential. Below threshold, the potential leads to bound states, above threshold it 
affects the cross section through final state interactions. It is, therefore, natural to take for the QCD 
case the coupling governing the QCD potential at the momentum scale involved in the rescattering. 
To relate ay to cr^g, we use 

a W s(Q 2 ) = a v (e+ 5 ^Q 2 ) 

In a similar way, BLM scale-fixing is adopted for the correction from hard gluon exchange. In the 
radiative correction, there still remain 0(a 2 s ) terms, identical to the radiative corrections for the 



1 + 2 ^ + 0(a v) 

IT 



(26) 
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Figure 2: Anisotropy for charmed, bottom and top quark production as a function of (3. Also shown 
is the Born prediction. We have assumed the effective quark masses m c = 1.7 GeV, nif, = 5 GeV 
and nit = 175 GeV. 



theory with a fixed coupling constant. With the same scheme convention for the coupling as above, 
one arrives at 

a v (4m 2 /3 2 /e) 



1 - e 



p ( 1 _ 4 «v(™V/ 6 r 
A = T (i 16 « v ( m 2 e 3/4) ^ -| _ e -x'„ av (4m 2 /3 2 ) 



where 



4ir a v (4m 2 /3 2 ) 
~3~ " 



, 4 7T ay (4m 2 /3 2 /e) 



(27) 



f28) 



/3 ' s 3 /3 
The anisotropy A is plotted in Fig. 2 versus the velocity /3 in the range 0.2 < (3 < 0.5 for charmed, 
bottom, and top quarks. For comparison, the tree level prediction is also shown. For charmed 
quarks, only (3 values above 0.4 are admitted in order to allow for the simultaneous production of 
DD and D*D. The charm prediction is particularly sensitive to the QCD parameters, since very 
low scales are accessible. Measurements of the anisotropy for e + e~ — ► cc thus have the potential of 
determining ay in the regime where perturbation theory begins to fail. 

The curves are based on an input value a^^~ 5 \M^) = 0.115. We use the two-loop beta-function 
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Figure 3: Sensitivity of the anisotropy A for (a) e + e —^bb and (b) e + e — ► cc to changes in 
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to evolve cnj^g to lower momenta and then used Eq. (26) to calculate ay(Q 2 ). To investigate the 
sensitivity of the predictions for bottom quarks, the input value for ajj^(M^) has been varied by 
±0.008 from the central value of 0.115. As demonstrated in Fig. 3a the variation of the anisotropy 
parameter amounts to about 10%, and could therefore be accessible experimentally. The charm 
predictions (see Fig. 3b) are even more sensitive. 

4. Conclusions 

An important consequence of heavy quark kinematics is that the production angle of a heavy hadron 
follows the direction of the parent heavy quark. This applies not only at Born approximation, but 
also after QCD corrections have been applied. In this paper we have shown that the anisotropy 
A{fl 2 ) in the cm angular distribution d<r(e + e — ► QQ)/dCl oc 1 + A cos 2 9 of heavy quarks produced 
near threshold is sensitive to the QCD coupling ay(Q 2 ) at specific scales determined by the quark 
relative momentum p cm = y/sfi. The coupling ay(Q 2 ) is the physical effective charge defined 
through heavy quark scattering. The predictions provide a connection between observables and 
thus are independent of theoretical conventions. 

An important feature of our analysis is the use of BLM scale-fixing, in which all higher-order 
corrections associated with the beta-function are resummed into the scale of the coupling. The 
resulting scale for ay(Q 2 ) corresponds to the mean gluon virtuality. In the case of the soft rescat- 
tering corrections to the 5-wave, the BLM scale is s (i 2 = p 2 m . One thus has sensitivity to the 
running coupling over a range of momentum transfers within the same experiment. The anisotropy 
measurement thus can provide a check on other determinations of ay(Q 2 ), e.g. from heavy quark 
lattice gauge theory, or from the conversion of cnj^g determinations to ay as given in Eq. (26). 

Our analysis also shows that the running coupling appears within the cross section with several 
different scales. This is particularly apparent at low (3 where the physical origin of the 0{a s ) 
corrections can be traced to gluons with different polarization and virtuality. 
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In principle, the anisotropy of r pairs produced in e + e~ — ► r + r~ could be used to measure the 
Pauli form factor F2(s) near threshold s ~ 4m^. A highly precise measurement of the anisotropy thus 
could provide a measurement of a fundamental parameter of the r lepton, its timelike anomalous 
magnetic moment. 
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